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1  ABSTRACT 
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3 

4  This  work  examines  the  applicability  of  a  previously-postulated  heuristic  model  for  the 

5  temporal  evolution  of  the  small-amplitude  tilt  of  a  tropical  cyclone-like  vortex  under 

6  vertical-shear  forcing  for  both  a  dry  and  cloudy  atmosphere.  The  heuristic  model  hinges 

7  on  the  existence  of  a  quasi-discrete  vortex -Rossby  wave  and  its  ability  to  represent  the 

8  coherent  precession  and  tilt-decay  of  a  stable  vortex  in  the  free-alignment  problem. 

9  Linearized  numerical  solutions  for  a  dry  and  cloudy  vortex  confirm  the  model  predictions 

10  that  an  increase  in  the  magnitude  of  the  radial  potential  vorticity  (PV)  gradient  within  the 

1 1  vortex  skirt  surrounding  the  core  yields  a  more  rapid  evolution  of  a  sheared  vortex 

12  towards  the  equilibrium,  left-of-shear  tilt  configuration.  However,  in  the  moist-neutral 

13  limit,  in  which  the  effective  static  stability  vanishes  in  rising  and  sinking  regions,  the 

14  heuristic  model  yields  a  poor  approximation  to  the  simulated  vortex  core  evolution,  but  a 

15  left-of-shear  tilt  of  the  near-core  vortex,  radially  beyond  the  heating  region,  remains  the 

16  preferred  long-time  solution.  Within  the  near-core  skirt  the  PV  perturbation  generated  by 

17  vertical  shearing  exhibits  continuous-spectrum  type  vortex-Rossby  waves,  features  that 

18  are  not  captured  by  the  heuristic  model.  Nevertheless,  the  heuristic  model  continues  to 

19  predict  the  rapid  vertical  alignment  and  equilibrium,  left-of-shear  tilt  configuration  of  the 

20  simulated  near-core  vortex  in  the  moist-neutral  limit. 
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1.  Introduction 


2  The  impact  of  environmental  vertical  wind  shear  on  tropical  cyclone  (TC) 

3  structure  and  intensity  depends  on  the  TC’s  dynamic  resilience  (Jones  1995;  Reasor  et  al. 

4  2004)  and  interaction  with  environmental  dry  air  (Simpson  and  Riehl  1958;  Kimball 

5  2006;  Riemer  and  Montgomery  2011),  both  of  which  influence  the  organization  of 

6  convection  and,  ultimately,  operation  of  the  TC  power  machine  (Emanuel  et  al.  2004; 

7  Riemer  et  al.  2010,  2013;  Tang  and  Emanuel  2010,  2012). 

8  The  problem  of  dynamic  resilience  focuses  on  the  ability  of  the  TC  to  maintain  a 

9  vertically-coherent  vortex  structure  as  it  experiences  vertical  shearing.  Jones  (1995) 

10  found  that  the  coupling  between  vertical  layers,  and  tendency  for  the  upper-  and  lower- 

1 1  level  potential  vorticity  (PV)  of  the  cyclonic  core  to  precess  into  an  upshear-tilted 

12  configuration,  restricts  the  development  of  vertical  tilt  that  would  otherwise  occur 

13  through  differential  advection.  For  small-amplitude  displacements  between  the  upper- 

14  and  lower-level  circulations,  Reasor  et  al.  (2004;  hereafter  RMG04)  developed  a  heuristic 

15  model  for  the  vertical  shear  forcing  of  vortex  tilt,  demonstrating  that  the  tilt  evolution 

16  depends  not  only  on  the  intrinsic  rotation  and  scales  of  the  flow  (e.g.,  Rossby  number  and 

17  Rossby  deformation  radius),  but  also  on  the  radial  distribution  of  PV  in  the  core  region  of 

18  the  vortex.  Here,  we  further  evaluate  this  heuristic  model. 

19  While  RMG04  showed  that  dramatically  different  vortex  evolution  could  be 

20  expected  for  different  vortex  profiles,  they  did  not  consider  a  set  of  controlled 

21  experiments  to  isolate  the  specific  impact  of  the  near-core  profile  on  the  sheared  vortex 

22  evolution.  To  understand  why  such  dependence  might  exist  in  the  first  place,  consider  the 

23  idealized  vortex  depicted  in  Fig.  1.  The  vortex  core  is  defined  as  the  innermost  region  of 
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the  vortex  bounded  approximately  by  the  radius  of  maximum  PV  gradient,  rc.  The  vortex 

2  skirt  immediately  surrounds  the  core  and  is  characterized  by  PV  values  much  less  than 

3  found  within  the  core.  We  will  alternately  refer  to  this  portion  of  the  vortex  as  the  near- 

4  core  region,  and  the  sum  of  the  core  and  skirt  as  the  core  region. 

5  When  a  vortex  is  tilted,  as  would  be  accomplished  by  a  burst  of  vertical  shear 

6  flow  from  the  environment,  the  PV  of  the  core  will  be  deformed.  In  the  limit  of  small  but 

7  finite-amplitude  tilt,  the  dynamics  of  the  core  can  be  meaningfully  described  by  the 

8  evolution  of  a  vortex -Rossby  wave  (VRW)  asymmetry,  referred  to  as  a  tilt  quasimode,  on 

9  a  background  axisymmetric  flow  (Reasor  and  Montgomery  2001).  Schecter  et  al.  (2002) 

10  alternatively  described  the  VRW  tilt  asymmetry  as  a  discrete  mode.  In  this  latter 

1 1  interpretation,  a  resonance  may  occur  at  the  critical  radius  where  the  precession 

12  frequency  of  the  tilt  mode  equals  the  angular  rotation  rate  of  the  background  flow.  As 

13  explained  in  Schecter  et  al.  (2002)  through  canonical  angular  momentum  (and 

14  equivalently,  angular  wave  activity)  conservation  arguments,  stirring  of  the  flow  near  the 

15  critical  radius  requires  a  damping  of  the  tilt  mode  if  the  radial  gradient  of  PV  at  rt  is 

16  negative.  Since  the  growth  in  PV  perturbation  amplitude  at  r*  depends  intrinsically  on 

17  how  sloped  the  background  PV  profile  is  there,  the  rate  at  which  this  damping  occurs, 

18  hence  the  vortex  resilience,  will  be  proportional  to  the  magnitude  of  the  PV  gradient 

19  within  the  vortex  skirt.  RMG04  extended  these  ideas  to  the  shear- forced  problem.  In 

20  section  3  we  revisit  their  heuristic  model  for  the  vortex  evolution  in  shear,  which  predicts 

21  a  similar  dependence  of  resilience  on  the  PV  gradient  at  r\. 

22  Since  moist  convection  pervades  tropical  cyclones,  especially  in  their  core  region, 

23  it  is  still  unclear  whether  the  resonant  damping  mechanism  using  dry  dynamics  governs 
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to  leading  order  the  vortex  resilience  in  shear.  For  one  thing,  the  static  stability  is  locally 

2  reduced  within  moist  convective  regions  of  the  flow  (e.g.,  Durran  and  Klemp  1982; 

3  Emanuel  et  al.  1987).  RMG04  observed  that  a  uniform  reduction  of  the  static  stability 

4  throughout  the  flow  can  lead  to  a  change  in  the  resilience  mechanism  from  one  governed 

5  by  resonant  damping  to  one  that  shows  greater  consistency  with  the  classic 

6  axisymmetrization  process  and  spiral  wind  up  (e.g.,  Melander  et  al.  1987;  Smith  and 

7  Montgomery  1995;  Montgomery  and  Kallenbach  1997).  They  did  not,  however, 

8  explicitly  consider  the  impacts  of  local  static  stability  reduction  on  the  resilience 

9  mechanism.  To  the  extent  that  the  resonant  damping  mechanism  may  still  govern  the 

10  resilience  of  a  moist-convective  vortex,  RMG04  reasoned  that  the  impact  of  heating 

1 1  would  be  to  reduce  the  local  deformation  radius,  thereby  increasing  the  azimuthal 

12  propagation  speed  of  the  tilt  mode.  An  increased  propagation  speed  would,  in  turn,  shift 

13  r*  inward  toward  the  core  where  the  PV  gradient  is  larger,  leading  to  more  effective 

14  damping  of  the  tilt  mode  (e.g.,  Schecter  and  Montgomery  2007;  see  their  Figs.  10  and  1 1 

15  and  accompanying  discussion).  In  a  recent  moist  numerical  study  of  TCs  in  shear,  Riemer 

16  et  al.  (2013)  documented  differences  in  vortex  tilt  evolution  potentially  related  to 

17  differences  in  the  initial  radial  vorticity  profile  outside  the  core.  To  investigate  this 

18  important  issue,  we  will  consider  the  sheared  vortex  evolution  under  a  simple 

19  parameterization  of  convective  heating  for  a  cloudy  vortex.  Here  “cloudy”  refers  to 

20  regions  of  non-precipitating,  saturated  air.  The  results  with  neutral  stratification  within  a 

21  limited  “eyewall”  region  will  show  that  the  heuristic  model  poorly  approximates  the 

22  overall  simulated  vortex  evolution,  calling  for  further  treatment  of  the  moist  problem  to 

23  fully  understand  the  mechanisms  for  TC  resilience. 
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1  The  paper  is  organized  as  follows.  Section  2  presents  the  numerical  model  with 

2  shear  forcing  and  parameterized  convective  heating.  The  heuristic  model  for  the  shear- 

3  forced  evolution  of  vortex  tilt  is  reviewed  in  section  3.  Idealized  numerical  simulations, 

4  which  address  the  relevance  of  the  near-core  profile,  are  presented  here.  Section  4  then 

5  considers  the  shear- forced  evolution  of  a  cloudy  vortex.  A  discussion  of  our  results  and 

6  conclusions  are  presented  in  section  5. 

7 

8  2.  Model 

9  Our  numerical  simulations  of  the  idealized  TC  in  shear  will  use  the  linear 

10  primitive  equation  model  described  in  RMG04,  modified  here  to  include  a 

1 1  parameterization  of  diabatic  heating.  The  model  assumes  a  stably  stratified  fluid  and 

12  employs  the  Boussinesq  approximation  on  an  /  plane.  The  model  is  formulated  in  a 

13  cylindrical  coordinate  system  (r,  cp,  z),  where  r  is  radius,  (p  is  azimuthal  angle,  and  z  is  the 

14  pseudoheight  vertical  coordinate  (Hoskins  and  Bretherton  1972),  and  is  bounded 

15  vertically  by  isothermal  rigid  lids  at  z  =  0  and  z  =  H.  Consistent  with  small  but  finite 

16  amplitude  departures  from  vertical  alignment,  the  governing  equations  are  linearized 

17  about  a  circular  (mean)  vortex  in  gradient  and  hydrostatic  balance.  The  mean  vortex  and 

18  quasi-static  flow  perturbation  (i.e.,  the  tilt  asymmetry)  are  then  decomposed  into 

19  barotropic  and  internal  baroclinic  modes.  As  in  RMG04,  for  simplicity,  we  consider  only 

20  initially  barotropic  mean  vortices.  Under  the  Boussinesq  approximation,  the  vertical 

21  structure  of  the  modes  comprising  the  geopotential  is  then  given  by  cos [\nnzjH],  where 

22  m  is  the  vertical  mode  number  and  H  is  the  physical  depth  of  the  vortex,  taken  to  be  10 

23  km  in  all  experiments. 
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1  The  tilt  asymmetry  is  Fourier  decomposed  in  azimuth.  The  prognostic  primitive 

2  equations  for  each  vertical  (in)  and  azimuthal  (n)  mode  subject  to  vertical  shear  forcing 

3  are  given  as  follows: 

4  f  4  + inQ  \mn  ('%  t)  -  &mn  +  =  -a(r)umn  +  F  ,  (1) 

\ot  J  or 

vmn  (r,  t )  +  rfumn  +  —  0mn  =  -a{r)vmn  +  F  ,  (2) 

r 

6  f  J;  +  in°  \k,n  (r,  t)+Y (r)cl  - -TT  (fUmn  )  +  =  ~cc{r%n  +  F ^  ,  (3) 

y ct  ^  [r  dr  r  \ 

7  where  umn ,  v„m ,  and  <pmn  are  the  radial  velocity,  tangential  velocity,  and  geopotential 

8  Fourier  amplitudes,  respectively.  The  basic  state  vortex  quantities,  £2  =v/r , 

9  %  =  f  +  2C2 ,  and  rf  =  f  +  C,  ,  are  the  azimuthal -mean  angular  velocity,  modified 

10  Coriolis  parameter,  and  absolute  vertical  vorticity,  respectively.  Here  C,  =  £ 1  d(rv)/dr  is 

11  the  relative  vertical  vorticity  and /is  the  constant  Coriolis  parameter,  taken  to  be  3.53  x 

12  10'5  s'1  (corresponding  to  14°  latitude)  in  all  experiments.  The  internal  gravity  wave  phase 

13  speed  cm  =  NH/mn ,  where  N  is  the  Brunt- Vaisala  frequency  and  is  assumed  to  take  a 

14  constant  value  of  1.22  x  10'“  s'  in  all  cases.  The  internal  gravity  wave  phase  speed  for  m 

15  =  1  is  approximately  40  m  s'1.  The  function  Y(r )  is  a  buoyancy  reduction  factor  (Schecter 

16  and  Montgomery  2007)  accounting  for  radial  variations  in  the  static  stability  in  the  core 

17  region,  and  is  further  motivated  below.  The  first  tenn  on  the  rhs  of  each  equation 

18  represents  Rayleigh  damping  within  a  “sponge”  ring  at  the  outennost  portion  of  the 

19  computational  domain  so  as  to  prevent  the  spurious  reflection  at  the  outer  boundary  of 

20  outward-propagating  inertia-gravity  waves  generated  from  the  vortex-shear  interaction  in 


8  .  yz 

- b  inQ 

dt 
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1  the  core  region.  The  sponge  ring  extends  inward  150  km  from  the  outer  boundary  with  a 

2  peak  damping  rate  of  1/300  s’1.  A  3000-km  outer  radius  domain  with  radial  grid  spacing 

3  Ar  =  1  km  is  used  in  the  experiments  of  section  3.  In  the  cloudy  vortex  experiments  of 

4  section  4,  the  domain  extends  to  1600  km  with  a  finer  grid  spacing  of  Ar  =  0.5  km. 

5  Additional  details  of  the  numerical  system  can  be  found  in  RMG04. 

6  The  mean  PV  for  a  barotropic  vortex,  q  =  N't]  ,  is  prescribed  and  the  associated 

7  tangential  velocity  field  is  assumed  in  gradient  balance: 


8 


(4) 


9  The  PV  perturbation  amplitude  derived  from  Eqs.  (l)-(3)  is  given  by 

(  \2 

,  „  A  ,  r2  r  W17T  — 2  ,C\ 

10  qmn=N  Cmn~  —  >1  </>,„„■  (5) 

V  H  J 


11  In  the  linear  framework  here,  a  tilted  vortex  is  approximated  by  adding  to  q  a  PV 


12  perturbation  with  m  =  n  =  1  and  radial  structure  proportional  to  the  radial  gradient  of  q  . 

13  This  approximation  to  a  tilt  asymmetry  takes  the  fonn  of  a  wavenumber- 1  dipole  with 

14  reflection  about  z  =  77/2.  Details  of  the  method  for  obtaining  quasi-balanced  perturbation 

15  asymmetric  wind  and  geopotential  fields  are  described  in  RMG04. 

16  The  second  term  on  the  rhs  of  Eqs.  (l)-(3)  is  a  linearized  forcing  term  that  arises 

17  with  the  presence  of  an  imposed  zonal  vertical  shear  flow,  which  is  assumed  weak 

18  compared  to  the  basic  state  swirling  flow.  For  simplicity,  the  zonal  shear  flow  is  chosen 

19  to  have  vertical  structure  of  the  first  internal  baroclinic  mode  so  that  only  vortex 

20  perturbations  with  m  =  n  =  1  are  forced.  In  cylindrical  coordinates,  the  radial  and 

2 1  azimuthal  velocity  components  of  the  shear  flow  forcing  are  expressed  as 
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7  To  represent  the  convective  heating/cooling  processes,  we  follow  the 

8  thermodynamically  reversible  formulation  of  Schecter  and  Montgomery  (2007)  and 

9  introduce  the  heating  into  the  equivalent  mass-continuity  Eq.  (3).  Under  the  assumption 

10  of  hydrostatic  motions,  they  derived  a  thermodynamic  equation  in  pressure  coordinates 

1 1  which  contains  expressions  for  the  moist  adiabatic  derivative  of  density  potential 

12  temperature  with  respect  to  pressure,  holding  moist  entropy  and  total  water  mixing  ratio 

13  constant,  in  saturated  and  unsaturated  air.  The  wave  equations  for  a  TC-like  vortex 

14  containing  both  clear  (i.e.,  unsaturated)  and  cloudy  (i.e.,  saturated)  regions  were  derived 

15  by  neglecting  the  product  of  the  vertical  motion  asymmetries  with  asymmetries  in  the 

16  thermodynamic  derivatives.  Although  mathematically  consistent,  the  accuracy  of  such  an 

17  approximation  in  the  context  of  the  tilt  asymmetry  evolution  has  not  yet  been  explored 

18  systematically.  In  the  resulting  thennodynamic  equation,  the  tenn  involving  the  static 

19  stability  is  modified  from  its  dry  counterpart  through  the  addition  of  a  buoyancy 

20  reduction  factor  7(r,  p),  which  accounts  for  the  azimuthally-averaged  clear  and  cloudy  air 

21  stability  properties  (Schecter  and  Montgomery  2007).  In  essence,  the  presence  of  cloudy 
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1  air  leads  to  a  reduction  of  the  azimuthal-mean  static  stability  (cf.  Patra  2004).  For  the 

2  barotropic  mean  vortices  considered  here,  we  assume  Y  is  independent  of  height.  Thus, 

3  the  thennodynamic  equation  of  RMG04  is  modified  by  replacing  the  constant  N  with 

4  Y(r)N  in  Eq.  (3).  In  the  case  of  a  clear  (non-cloudy)  vortex,  Y  — >  1  and  the  dry 

5  thermodynamic  equation  of  RMG04  is  recovered.  In  the  eyewall  of  a  hurricane,  we  can 

6  anticipate  values  of  Y  approaching  zero  (Schecter  and  Montgomery  2007;  see  their  Fig. 

7  2). 

8  For  the  cases  examined  here  (Table  1),  simulation  experiments  are  performed 

9  with  and  without  vertical  wind  shear.  In  the  absence  of  shear,  the  vortex  is  given  an  initial 

10  tilt  as  described  above.  The  ensuing  vortex  evolution  is  referred  to  here  as  ‘free 

11  alignment.’  In  the  experiments  of  section  3,  a  dry  Rankine-with-skirt  (RWS)  vortex  is 

12  used.  The  vortex  radial  profiles  shown  in  Fig.  2a  are  similar  except  for  different  radial 

13  gradients  of  PV  within  the  skirt  immediately  outside  the  core.  The  vortex  of  the  control 

14  experiment  RWS1  exhibits  the  weakest  gradient  within  the  skirt.  Compared  to  RWS1,  the 

15  vortices  of  RWS2  and  RWS3  have  radial  PV  gradients  that  are  33%  and  100%  greater, 

16  respectively.  The  RWS2  and  RWS3  vortices  do  have  greater  circulation  outside  the  core, 

17  but  in  all  cases  the  contribution  of  the  skirt  to  the  total  circulation  in  the  core  region  is 

18  small.  Inside  r  =  69  km  (the  critical  radius  estimated  below  in  section  3b)  of  the  RWS1, 

19  RWS2,  and  RWS3  vortices,  the  skirt  contributes  approximately  3.8%,  4.4%,  and  5.7%, 

20  respectively,  to  the  total  circulation. 

21  The  vortex  radial  profile  shown  in  Fig.  2b  is  used  in  the  “cloudy  vortex” 

22  experiments  of  section  4.  It  is  based  upon  low-level  radar  and  fight-level  observations  of 

23  Hurricane  Guillermo  (1997;  Reasor  and  Eastin  2012).  A  monotonic  approximation  to  the 
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1  observed  PV  profile  is  used  here  to  exclude  the  possibility  of  ‘dry’  exponential  barotropic 

2  instability  (Schubert  et  al.  1999)  and  the  trochoidal  wobble-mode  instability  (Nolan  et  al. 

3  200 1)1.  In  a  series  of  experiments,  discussed  further  in  section  4,  the  static  stability 

4  within  the  region  of  strongest  winds  is  varied  from  the  limit  of  no  heating  (OBS1)  to  the 

5  limit  of  neutral  values  (OBS5). 

6 

7  3.  Dependence  of  Vortex  Evolution  on  Near-core  PV  Profile 

8  a.  A  heuristic  model  revisited 

9  For  vertically  sheared  vortices  like  those  described  in  section  2,  the  phase  and 

10  amplitude  evolution  of  the  tilt  asymmetry  may  be  understood  through  a  simple  heuristic 

1 1  model.  As  in  RMG04,  we  assume  that  the  vortex  admits  a  free-alignment  solution  in 

12  which  the  tilt  asymmetry  takes  the  form  of  a  damped  discrete  VRW  mode  (Schecter  and 

13  Montgomery  2003).  The  PV  perturbation  may  then  be  expressed  (for  m  =  n  =  1)  as 

A|  +  c.c„  yd<  0,  (8) 

15  where  the  complex  amplitude  A(r )  provides  the  radial  structure  of  the  tilt  asymmetry,  ycl 

16  is  the  exponential  (Landau)  damping  rate  (Briggs  et  al.  1970),  cop  is  the  vortex 

17  precession  frequency  and  “c.c.”  denotes  complex  conjugate.  Here  yd  is  defined  as 

18  negative  for  the  decaying  mode.  The  temporal  part  of  (8)  is  the  solution  to  a  single- 

19  component  damped  harmonic  oscillator  equation  of  the  fonn 

20  ^-  +  (ico p-ycl)q'  =  0.  (9) 

dt 

1  Exponential  growth  of  perturbations  may  still  occur  through  the  excitation  of  inertia-gravity  waves  by  the 
rotating  tilt  asymmetry  (Schecter  and  Montgomery  2004,  2006). 


q'(r,<p,z,t)  =  A(r)e 


Td*  J\V>-(0Pt)  , 
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1  Equation  (9)  may  be  motivated  physically  from  a  re-casting  of  the  linearized  PV  equation 

2  without  shear  forcing: 

3  — +  Q-^-  q'(r,(p,z,t)+ u' —  =  0.  (10) 

[dt  dq>)  K  )  dr 

4  Note  that  Eq.  (10)  does  not  explicitly  contain  the  PV  source  tenn  derived  from  the 

5  diabatic  heating  in  the  system  (l)-(3).  In  the  heuristic  model  given  by  Eq.  (9)  the 

6  response  to  heating  is  represented  through  a  modification  of  the  mode’s  propagation.  The 

7  PV  equation  (10)  is  consistent  with  this  approximation.  The  heuristic  model  for  the  shear- 

8  forced  system  is  derived  by  first  introducing  to  the  rhs  of  Eq.  (10)  the  linearized  shear 

9  forcing  consistent  with  (7), 


10  Fsr  =  -U i~  +  “y- ~  >cos(^)cosf^l,  m  =  1.  (11) 

[dr  lrG  f)  VH ) 

1 1  Here  l  r  G  =  NH/k  f  is  the  “global”  Rossby  defonnation  radius  of  the  flow  (for  m  =  1). 


12  The  forced,  damped  core-mode  model  derived  below  is  partly  justified  by  the  strong 

13  projection  of  the  above  shear  forcing  onto  the  tilt  mode  wavefunction  (cf.  Schecter  et  al. 

14  2000).  The  solution  to  Eq.  (9)  with  the  forcing  (11)  added  to  the  rhs  is 


q'{r,q>,z,i)  =  A{r)eYc,t ^  cos  —  - 

V  H  ) 


u 

r 

dq 

iap-yc, 

dr 

< 

r —  —  —  l(D 

q  v  q  e r 


16  The  second  term  in  the  shear  forcing  (11),  and  which  appears  in  the  solution  (12),  is 

17  associated  with  the  horizontal  advection  of  vertical  shear-flow  PV  by  the  basic  state 

18  vortex  flow.  Within  the  core  region  of  the  vortices  in  Fig.  2,  a  scale  analysis  and 

19  calculation  show  that  this  tenn  is  at  least  an  order  of  magnitude  smaller  than  the  first  term 

20  representing  the  differential  advection  of  basic  state  vortex  PV  by  the  vertical  shear  flow. 
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1  Thus,  for  the  vortices  here,  the  solution  may  be  further  simplified  by  neglecting  this 

2  second  tenn  of  the  shear  forcing: 


4  Note  an  implicit  assumption  in  the  derivation  of  (12)  and  (13)  is  that  ycl  and  cop  of  the 

5  free-alignment  solution  (8)  continue  to  govern  the  forced  evolution  of  the  vortex.  In  the 

6  fully  nonlinear  problem  this  should  be  a  valid  approximation  provided  the  shear  forcing 

7  is  relatively  weak. 

8  Assuming  the  initial  vortex  is  vertically  aligned  and  «  1  (Schecter  and 

9  Montgomery  2003),  it  follows  analytically  from  (13)  that  in  the  large-time  limit  the 
10  forced  solution  will  be  governed  approximately  by 


q'(r,<p,z,t )«  U  ^  eTc,t : 
cop  dr 


f)-sin^  c 


12  This  heuristic  model  predicts  that  small  but  finite  damping  results  in  a  long-time  solution 

13  of  fixed  tilt  as  measured  by  the  amplitude  and  phase  of  the  PV  asymmetry.  Since  dq/dr 

14  is  negative  outside  the  core  for  “real-world”  TC  vortices  (Mallen  et  al.  2005),  the 

15  equilibrium  vortex  tilt  is  directed  90°  left  of  shear.  Furthermore,  since  ycI  is  proportional 


16  to  the  PV  gradient  at  the  critical  radius  r*  of  the  discrete  VRW  (from  the  free-alignment 


17  solution),  the  resilience  of  the  shear-forced  vortex  should  depend  upon  the  near-core 

18  tangential  wind  profile.  We  validate  this  prediction  of  the  heuristic  model  through 

19  explicit  numerical  simulations  discussed  below.  The  heuristic  model  solution  (14)  in  the 

20  limit  of  ycl=  0  has  recently  been  confirmed  by  Paschke  et  al.  (2012)  for  the  fully 
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nonlinear  problem  through  the  method  of  matched  asymptotic  expansions.  Note  that  in 
this  limit  (14)  may  be  re-expressed  as 


2 

3 
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2U  dq  . 

- sin| 

(op  dr 


f  (Opt  ^ 
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cos 
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(Opt 


cos 


\H  j 


(15) 


Thus,  according  to  the  heuristic  model  in  the  absence  of  damping,  the  shear-forced  tilt 
asymmetry  should  precess  with  half  the  frequency  of  the  intrinsic  tilt  mode  in  the  free- 
alignment  case. 


b.  Heuristic  model  validation 

As  noted  in  section  2,  the  RWS  vortices  used  in  RWS1-3  differ  primarily  in  the 
gradient  of  PV  within  the  near-core  skirt  (Fig.  2a).  The  intrinsic  resilience  of  each  vortex 
is  explored  first  by  adding  to  the  symmetric  vortex  the  balanced  PV  perturbation 
approximating  the  m  =  n  =  1  tilt  asymmetry.  For  brevity,  we  will  denote  the  geopotential 
Fourier  amplitude  for  m  =  n  =  1  by  |0X|  and  the  corresponding  three-dimensional 
perturbation  value  by  <p1.  A  similar  convention  will  be  employed  for  all  perturbation 
variables.  A  time  series  of  the  geopotential  perturbation  amplitude  evaluated  at  the  core 
radius  rc  is  shown  in  Fig.  3.  Here  and  elsewhere,  time  is  expressed  in  units  of  circulation 
period,  defined  here  as  Inj Q(rc ) . 

From  Fig.  3,  during  the  first  circulation  period,  the  asymmetry  adjusts  as  it  takes 
the  fonn  of  the  tilt  mode.  Thereafter,  the  decay  of  the  tilt  asymmetry  is  approximately 
exponential"  for  multiple  circulation  periods,  consistent  with  prior  theoretical  analyses 
(Schecter  et  al.  2002;  Schecter  and  Montgomery  2003).  While  semi-analytical  methods 

2  The  low-frequency  oscillations  in  amplitude  result  from  the  initialization  of  PV  perturbation  in  the  critical 
layer  of  the  tilt  mode  (Schecter  and  Montgomery  2006).  The  oscillations  are  thus  not  a  fundamental  feature 
of  the  solution. 
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exist  to  compute  cop  and  ycl  (e.g.,  Schecter  et  al.  2002,  Appendix),  we  estimate  the 

2  values  here  by  constructing  linear  best  fits  to  the  simulated  geopotential  phase  and 

3  (logarithmically-scaled)  amplitude  following  the  initial  adjustment  period.  It  is  important 

4  to  recognize,  however,  that  the  numerical  solution  reflects  both  the  effects  of  resonant 

5  damping  and  tilt  asymmetry  growth  through  inertia-gravity  wave  (IGW)  pumping 

6  (Schecter  and  Montgomery  2004,  2006).  Thus,  the  decay  rate  estimated  here  is  actually 

7  the  sum  of  the  (negative)  resonant  damping  rate  ycl  and  (positive)  exponential  IGW 

8  pumping  rate  yrad  .  We  define  y  =  yrad  +  ycl  as  the  net  decay  rate  to  avoid  confusion. 

9  Table  2  lists  the  estimated  values  of  r„  cop ,  and  y  for  each  of  the  three  RWS 

10  experiments.  The  value  of  cop  increases  slightly  from  RWS1  to  RWS3  due  to  the 

1 1  increased  amount  of  circulation  immediately  outside  the  core.  This  increase  has 

12  negligible  impact  on  the  value  of  r„.  The  value  of  y ,  however,  increases  by  an  order  of 

13  magnitude  with  the  doubling  of  the  PV  gradient  within  the  skirt.  Schecter  and 

14  Montgomery  (2003)  showed  that  the  resonant  damping  rate  ycl  is  directly  proportional  to 

15  the  PV  gradient  at  r„.  In  the  absence  of  IGW  pumping,  then,  the  decay  rate  ought  to 

16  approximately  double  from  RWS1  to  RWS3.  Note  that  if  the  IGW  pumping  were  to 

17  largely  offset  the  resonant  damping  in  RWS1,  but  not  in  RWS3,  then  y  could  in  principle 

18  increase  by  a  factor  substantially  larger  than  2.  The  weak  decay  observed  in  Fig.  3  for 

19  RWS1  suggests  that  this  is  a  plausible  explanation  for  the  order-of-magnitude  increase  in 

20  net  decay  rate  . 


3  Estimates  of  /rad  and  Ycl  using  the  simulated  solution  amplitudes  in  the  analytical  expressions  of 
Schecter  and  Montgomery  (2004,  Eqs.  61  and  68)  suggest  that  yd  is  only  ~7.5%  larger  than  yrad  in 
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1  The  same  vortices  with  no  initial  vertical  tilt  are  then  subjected  to  unidirectional 

2  shear.  In  all  cases  we  take  U  =  -2  m  s’1,  yielding  4  m  s’1  of  deep-layer  westerly  shear. 

3  Figure  4  shows  the  PV  perturbation  amplitude  measured  at  rc  and  the  direction  of  tilt, 

4  defined  here  as  the  PV  perturbation  crest  azimuth  at  z  =  H,  for  each  case4.  In  RWS1  the 

5  vortex  tilts,  precesses  upshear,  and  realigns  on  a  time  scale  of  5  circulation  periods.  This 

6  evolution  is  then  repeated.  The  direction  and  amplitude  show  signs  of  a  slow  approach  to 

7  a  quasi-equilibrium  state.  The  rotation  frequency  estimated  from  the  PV  asymmetry 

8  phase  in  Fig.  4a  is  approximately  half  the  value  of  (O  for  this  weakly-damped  case,  in 

9  agreement  with  the  heuristic  model  prediction  given  by  (15).  The  increased  skirt  PV 

10  gradient  in  RWS2  and  RWS3  does  in  fact  lead  to  a  more  rapid  evolution  of  the  vortex  tilt 

1 1  towards  a  quasi-equilibrium  left-of-shear  orientation.  After  the  initial  downshear  tilting, 

12  the  vortex  in  RWS3  oscillates  about  the  left-of-shear  direction,  never  departing  from  this 

13  orientation  by  more  than  45°. 

14  The  PV  amplitude  and  tilt  direction  from  the  heuristic  model  solution  (14),  using 

15  the  values  of  cop  and  y  from  Table  2,  are  also  shown  in  Fig.  4.  The  net  decay  rate  is  used 

16  in  place  of  the  resonant  damping  rate  in  (14)  since  it  is  the  simulated  free-alignment 

17  decay  that  is  relevant  here.  The  peak  amplitude  value  predicted  by  (14)  should  be  most 

18  accurate  (i.e.,  reproduce  the  simulated  value)  when  the  vortex  closely  approximates  a 

19  single-core -plus-weak-skirt  structure.  While  we  have  verified  that  this  is  indeed  the  case 

20  for  the  RWS  vortices  here  (not  shown),  non-negligible  quantitative  differences  in  peak 


RWS1  (not  shown).  In  RWS3  the  estimated  value  of  yc/  is  more  than  double  yrad  so  that  the  order-of- 
magnitude  increase  in  y  from  RWS1  to  RWS3  is  largely  accounted  for. 

4  Tilt  evolution  is  depicted  using  PV  in  the  forced  experiments  to  parallel  the  heuristic  model  development 
in  section  3a.  Geopotential,  which  lacks  the  higher  frequency  oscillations  evident  in  the  PV  field,  is  better 
suited  for  estimating  the  decay  rate  and  precession  frequency  in  the  free-alignment  experiments. 
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1  amplitude  for  the  broader  vortex  of  section  4  are  found.  Thus,  for  consistency,  the  peak 

2  amplitude  from  the  heuristic  model  is  scaled  so  as  to  match  the  initial  peak  simulated 

3  value  for  all  experiments.  Despite  being  based  upon  the  approximate  solution  (8)  using 

4  the  vortex  free-alignment  response,  the  heuristic  model  describes  extraordinarily  well  the 

5  simulated  evolution  of  the  tilt  asymmetry.  Effects  related  to  the  horizontal  advection  of 

6  the  vertical  shear-flow  PV  by  the  basic  state  vortex  flow  (i.e.,  the  second  tenn  on  the  rhs 

7  of  (1 1))  are,  as  discussed  above,  subdominant  to  the  impacts  of  differential  advection  of 

8  basic  state  vortex  PV  by  the  vertical  shear  flow.  At  long  times,  the  radial  structure  of  the 

9  simulated  PV  perturbation  is  proportional  to  dq/dr  (not  shown),  as  anticipated  by  the 

10  heuristic  model  solution  (14).  Therefore,  it  follows  that  the  enhanced  resilience  of  the 

11  RWS2  and  RWS3  sheared  vortices  can  be  attributed  directly  to  the  increased  radial 

12  gradient  of  PV  within  the  skirt  and  its  impact  on  the  resonant  damping  of  the  tilt  mode. 

13  When  one  considers  the  nature  of  the  resonant  damping  mechanism,  that  damping 

14  occurs  due  to  the  resonance  of  a  propagating  core  VRW  with  the  surrounding  flow  at  the 

15  critical  radius  of  the  core  mode,  it  might  at  first  sight  seem  surprising  that  the  forced 

16  solution,  which  exhibits  a  time-dependant  rate  of  phase  propagation  (Fig.  4),  should  be 

17  closely  reproduced  by  a  heuristic  model  invoking  values  of  cop  and  y  based  upon  a 

18  constantly-propagating  mode. 

19  To  gain  a  more  complete  understanding  of  the  forced  solution,  we  begin  by 

20  examining  the  development  of  the  linear  critical  layer.  Figure  5  a  shows  the  azimuthal 

21  structure  of  the  PV  perturbation  within  a  radial  band  centered  on  the  estimated  critical 

22  radius  for  the  RWS3  free-alignment  experiment.  As  basic  state  PV  within  the  skirt  is 

23  advected  radially,  the  PV  perturbation  grows  within  the  critical  layer  and  perturbation 
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1  growth  here  is  most  effective  near  the  location  of  the  resonance.  Conservation  of  global 

2  angular  wave  activity  requires  that  the  amplitude  of  the  tilt  mode  must  then  decay  (e.g., 

3  Schecter  et  al.  2002). 

4  A  similar  illustration  of  the  flow  evolution  is  shown  in  Fig.  5b  for  the  shear- 

5  forced  RWS3  experiment.  A  critical  layer  appears  to  develop  near  the  same  radial 

6  location  observed  in  the  free-alignment  experiment,  suggesting  that  the  resonant 

7  interaction  of  a  constantly-propagating  mode  continues  to  play  a  role  in  the  linear  tilt 

8  evolution  in  the  forced  problem.  The  temporal  evolution  of  the  PV  perturbation  within 

9  the  skirt  in  the  shear-forced  RWS3  experiment  is  shown  in  greater  detail  in  Fig.  6.  Over 

10  the  duration  of  the  simulation,  the  peak  amplitude  of  the  PV  perturbation  within  the  skirt 

1 1  grows.  The  corresponding  decrease  in  the  peak  amplitude  of  the  core  tilt  asymmetry  (Fig. 

12  4c)  agrees  qualitatively  with  the  predictions  of  the  resonant  damping  theory.  The 

13  amplitude  of  the  PV  perturbation  within  the  skirt,  however,  is  modulated  by  the  shear 

14  forcing  just  as  it  is  in  the  core  (Fig.  4c). 

15  If  the  critical  layer  stirring  were  instead  tied  to  the  instantaneous  rotation 

16  frequency  of  the  shear-forced  solution,  one  would  expect  to  observe  significant  radial 

17  variation  in  the  location  of  rt  with  time.  Except  for  a  brief  interval  of  time  near  5 

1 8  circulation  periods,  the  estimated  critical  radius  using  the  rotation  frequency  of  the  shear- 

19  forced  solution  lies  outside  90-km  radius  at  all  times  during  the  RWS3  simulation  (not 

20  shown).  The  observed  peak  in  PV  perturbation  amplitude  within  the  skirt,  by  contrast,  is 

21  consistently  located  near  the  value  of  r,.  from  the  free-alignment  simulation.  Thus,  we 

22  conclude  that  resonant  damping  in  the  linear  shear-forced  problem  continues  to  be 

23  governed  essentially  by  the  intrinsic  tilt  mode  identified  in  the  free  alignment  experiment. 
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1  The  heuristic  model  predicts  that  the  amplitude  of  this  rotating  component  eventually 

2  decays  to  zero,  at  which  point  resonant  damping  would  cease.  In  this  limit,  the  vortex  has 

3  achieved  an  equilibrium  tilt  90°  left  of  shear. 

4 

5  4.  “Cloudy  Vortex”  Experiments 

6  When  a  TC  encounters  vertical  wind  shear,  observations  confirm  that  deep 

7  convection  within  the  near-core  region  organizes  into  a  distinctly  asymmetric  pattern 

8  with  maximum  ascent,  precipitation,  and  lightning  frequency  downshear  left  (e.g.,  Reasor 

9  et  al.  2000;  Corbosiero  and  Molinari  2003;  Reasor  et  al.  2013).  Several  plausible 

10  explanations  have  been  put  forth  in  the  literature  for  how  this  occurs,  but  yet  there  is  no 

1 1  clear  consensus  as  to  which  mechanism  is  predominant.  One  set  of  theories  argues  that 

12  the  balanced,  adiabatic  dynamics  of  the  vortex-shear  interaction  determines  regions 

13  favored  for  ascent.  Diabatic  heating  is  then  organized  within  the  envelope  of  upward 

14  mesoscale  ascent  in  both  conditionally  stable  and  unstable  environments  (e.g.,  Davis  et 

15  al.  2008;  Riemer  et  al.  2010).  The  parameterization  of  heating  in  the  system  of  equations 

16  (l)-(3)  follows  this  conceptual  view  by  assuming  that  the  diabatic  heating  (cooling)  is 

17  coupled  to  the  asymmetric  wavenumber- 1  ascent  (descent)  associated  with  the  dry  vortex 

18  dynamics.  The  net  effect  of  this  coupling  is  to  locally  reduce  the  azimuthal-mean  static 

19  stability  within  a  radial  band.  A  similar  approach  was  taken  by  Patra  (2004)  in  their 

20  examination  of  the  convective  heating’s  impact  on  the  resilience  of  a  sheared  TC  vortex. 

21  Figure  7  shows  the  radial  profiles  of  buoyancy  reduction  factor  Y  for  the  cloudy 

22  vortex  experiments  performed  here.  To  establish  the  basic  resilience  properties  of  the 

23  observationally-based  vortex,  the  first  experiment,  OBS1,  involves  no  heating  whatsoever 
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(i.e.,  7=1  everywhere).  The  remaining  experiments  specify  7  <  1  at  the  outer  edge  of  the 

2  core,  rc  defined  herein  by  the  radius  of  maximum  PV  gradient.  (This  specification  is 

3  consistent  with  the  finding  of  Schecter  and  Montgomery  (2007)  using  a  non-hydrostatic 

4  mesoscale  model  that  the  minimum  7  coincides  with  the  maximum  azimuthal-mean 

5  tangential  wind.)  In  OBS2-4,  the  value  of  7  is  progressively  reduced  locally  and 

6  asymptotes  to  the  OBS1  value  within  the  “eye”  and  vortex  skirt.  Experiment  OBS5  is 

7  similar,  except  7  is  locally  reduced  to  zero. 

8  The  free  alignment  of  the  vortex  in  each  case,  measured  at  rc,  is  shown  in  Fig.  8a. 

9  Following  the  initial  adjustment  of  the  perturbation’s  radial  structure  in  OBS1,  as  in  the 

10  experiments  of  section  3,  a  clear  exponential  decay  of  the  tilt  asymmetry  is  observed.  The 

1 1  critical  radius  r*  ~  3 rc  in  this  case  (Table  2).  The  local  reduction  of  static  stability  in 

12  OBS2  and  OBS3  results  in  an  increase  in  the  value  of  (O  and  corresponding  inward  shift 

13  of  r„.  Also,  y  in  OBS2  (OBS3)  is  more  than  a  factor  of  two  (three)  larger  than  in  OBS1. 

14  A  further  reduction  of  the  local  static  stability  to  near-neutral  values  in  OBS4  results  in  a 

15  tilt  asymmetry  decay  that  is  not  clearly  exponential.  When  the  tilt  asymmetry  is  measured 

16  instead  at  r  =  40  km,  radially  outside  the  region  of  peak  heating,  exponential  decay  is 

17  observed  in  OBS4  following  the  initial  adjustment  (Fig.  8b)  at  a  rate  five  times  larger 

18  than  in  OBS1.  The  different  solution  behavior  inside  and  outside  the  heating  region  for 

19  large  reductions  of  static  stability  suggests  that  tilt  defined  via  conventional  methods  like 

20  PV  centroid  will  have  a  strong  radial  scale  dependence  in  the  TC.  In  OBS5,  where  the 

21  static  stability  is  reduced  locally  to  zero,  exponential  decay  of  the  tilt  asymmetry  is  not 

22  readily  observed  within  or  radially  outside  the  region  of  heating.  Without  employing  the 

23  previously-noted  semi-analytical  methods  for  computing  the  precession  frequency  and 
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damping  rate,  the  contribution  of  the  tilt  mode  to  the  tilt  asymmetry  evolution  cannot  be 

2  clarified  in  this  limiting  case.  Here  it  is  not  essential  that  the  weak  excitation  of  the  tilt 

3  mode  be  fully  quantified.  We  refer  the  reader  to  a  recent  study  by  Schecter  (2015)  which 

4  examines  the  issue  of  tilt  mode  excitation  in  the  context  of  convection-permitting  TC 

5  simulations. 

6  The  observation-based  vortex  is  next  subjected  to  unidirectional  shear  with  U  =  - 

7  2.75  m  s'1,  yielding  5.5  m  s'1  of  deep-layer  westerly  shear.  Figure  9  shows  the  PV 

8  perturbation  amplitude  measured  at  rc  (with  the  exception  of  OBS4-5,  where  it  is 

9  measured  at  r  =  40  km)  and  the  direction  of  tilt  in  each  case.  For  reference,  1  circulation 

10  period  is  equivalent  to  ~45  min.  Within  a  few  circulation  periods  the  tilt  asymmetry  in 

11  OBS1  (Fig.  9a)  approaches  a  quasi-equilibrium  configuration  in  which  the  tilt  direction 

12  varies  by  no  more  than  45°  from  the  left-of-shear  direction  and  the  tilt  amplitude 

13  gradually  asymptotes  to  a  steady  value.  This  long-term  resistance  to  shearing  reinforces 

14  the  idea  promoted  by  RMG04  and  Jones  (2004)  that  mechanisms  intrinsic  to  the  adiabatic 

15  dynamics  are  capable  of  maintaining  the  vertical  alignment  of  a  TC-like  vortex  on 

16  meteorologically-relevant  timescales5.  Furthermore,  the  close  agreement  between  the 

17  simulated  and  heuristic  model  solutions  supports  the  relevance  of  the  resonant  damping 

18  mechanism  to  vortices  more  closely  resembling  TCs.  The  simulated  solution  shows 

19  higher- frequency  hannonics  and  a  growing  departure  from  the  heuristic  model  solution 

20  with  time.  These  hannonics  are  not  evident  in  the  RWS  cases  of  section  3,  indicating  that 


5  Whether  the  vortex  is  able  to  survive  the  initial  tilting  in  the  fully  nonlinear  problem  is  not  formally 
addressed  here.  RMG04  argued  that  a  necessary  condition  required  for  the  vortex  to  resist  the  shearing  is 

revealed  by  comparing  the  free-alignment  value  of  cop  to  the  differential  advection  rate,  t~  ~U/L. 
According  to  this  perspective,  a  vortex  will  tend  to  remain  vertically  coherent  when  the  intrinsic  precession 
frequency  exceeds  r”1  .  For  the  OBS  experiments,  r^1  ~  10'4  s'1.  Since  cop  is  at  least  a  factor  of  5  larger  in 
each  of  the  OBS  cases,  we  may  anticipate  a  nonlinearly  resilient  vortex. 
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1  some  assumptions  and  approximations  of  the  single-core-plus-weak-skirt  heuristic  model 

2  start  to  become  tenuous  for  broader  radial  distributions  of  PV. 

3  The  shear-forced  tilt  asymmetries  in  OBS2  (Fig.  9b)  and  OBS3  (Fig.  9c)  with 

4  moderate  local  reduction  of  Y  behave  similarly  to  that  in  the  no-heating  experiment 

5  OBS1.  As  anticipated  from  the  discussion  summarized  in  the  Introduction,  the  resonant 

6  damping  is  enhanced  through  the  local  reduction  in  Rossby  deformation  radius  and 

7  inward  shifting  of  the  critical  radius  to  a  region  of  greater  PV  gradient.  As  a  consequence, 

8  the  vortex  achieves  the  quasi-equilibrium  tilt  configuration  over  a  shorter  time  period  and 

9  with  smaller  final  tilt  magnitude  than  without  heating. 

10  In  OBS4  (Fig.  9d)  and  OBS5  (Fig.  9e)  the  shear- forced  vortex  evolution  for  the 

1 1  transition  to  locally  neutral  static  stability  values  in  the  core  is  examined.  Because  the 

12  free-alignment  solution  shows  little  evidence  of  damped-mode  behavior  within  the  region 

13  of  peak  heating  in  these  cases,  the  tilt  magnitude  and  direction  are  instead  measured 

14  outside  the  region  of  heating  at  the  outer  edge  of  the  core.  In  contrast  to  the  evolution  in 

15  OBS1-3,  the  tilt  magnitude  in  OBS4  almost  immediately  exhibits  high-frequency, 

16  constant-amplitude  oscillations  about  a  quasi-equilibrium  value.  The  tilt  direction 

17  behaves  similarly.  Note,  however,  that  the  longer-timescale  variations  in  magnitude  and 

18  direction  predicted  by  the  heuristic  model  over  the  first  10  circulation  periods  are  still 

19  reflected  in  the  simulated  solution.  The  differences  between  the  simulated  tilt  evolution  in 

20  OBS4  and  OBS5  are  subtle.  Due  to  our  inability  to  theoretically  estimate  cop  and  y  in  the 

21  OBS5  free-alignment  case,  we  simply  reproduce  the  heuristic  model  solution  from  OBS4 

22  in  Fig.  9e  as  an  approximation.  As  in  OBS4,  the  simulated  solution  continues  to  reflect 

23  the  longer-time  oscillations  in  tilt  magnitude  and  direction  from  the  heuristic  model. 
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1  Further  insight  into  the  impact  of  heating  on  the  shear-forced  tilt  is  obtained  by 

2  comparing  the  PV  and  vertical  velocity  perturbation  structure  from  OBS1  (Fig.  10)  with 

3  that  from  OBS5  (Fig.  11)  at  select  times  during  the  early  evolution.  Only  the  ascending 

4  contribution  to  the  wavenumber- 1  vertical  motion  asymmetry  is  shown  for  clarity. 

5  Differential  advection  of  the  symmetric  vortex  by  the  vertical  shear  flow  at  early  times 

6  produces  a  broad-scale  PV  dipole  in  the  direction  of  the  imposed  shear  vector  (Fig.  10a 

7  and  11a).  According  to  balance  arguments  of  Jones  (1995),  ascent  (descent)  is  initially 

8  produced  downshear  (upshear)  as  the  vortex  attempts  to  maintain  thermal  wind  balance 

9  on  account  of  the  tilted  wind  field.  Since  the  parameterized  diabatic  heating  is 

10  proportional  to  the  vertical  motion,  it  can  be  readily  shown  using  the  definition  of  PV 

1 1  perturbation  for  the  Boussinesq  system  linearized  about  a  barotropic  mean  vortex, 

12  q’  =  N^'  +  Wt£,  (16) 

cz 


17  Thus,  within  the  “eyewall”  region  of  lowest  Y  in  OBS5,  positive  (negative)  PV 

18  perturbation  is  generated  downshear  (upshear)  at  low  levels  as  can  be  inferred  from  Fig. 

19  11a.  The  diabatic  generation  of  PV  perturbation  occurs  within  the  system-scale  tilt 

20  asymmetry  produced  by  differential  advection,  effectively  splitting  the  initial  tilt 

2 1  asymmetry. 


23 


1 


The  subsequent  vortex  tilt  evolution  is  readily  discerned  in  the  dry  limit 

2  configuration  of  OB  SI  given  the  discrete  structure  of  the  PV  asymmetry  throughout  the 

3  domain.  As  the  tilt  rotates  cyclonically  and  increases  in  magnitude,  the  peak  ascent  in 

4  OBS 1  shifts  to  the  right-of-tilt  direction  and  maintains  this  orientation  as  the  tilt  oscillates 

5  about  the  left-of-shear  direction  (Fig.  lOb-d).  Jones  (1995)  described  this  phase 

6  relationship  between  the  vertical  motion  asymmetry  and  tilt  as  the  balanced-dynamical 

7  response  of  vortex  flow  moving  over  the  distorted  isentropes  of  the  tilted  vortex.  The 

8  significant  local  reduction  of  azimuthal-mean  static  stability  in  OBS5  yields  a  more 

9  complex  PV  and  vertical  motion  structure  within  the  core  region  following  the  initial 

10  vortex  tilting  (Fig.  llb-d).  Inside  the  moist-neutral  “eyewall”  region  from  24-  to  28-km 

1 1  radius,  the  PV  perturbation  is  comprised  of  fine-scale  banded  structures.  Animations  of 

12  these  features  indicate  that  they  retrogress,  propagating  clockwise  (not  shown).  The 

13  vortex  tilt  on  this  sub-vortex  scale  is  difficult  to  assess  in  OBS5  without  employing  some 

14  quantitative  method  (e.g.,  PV  centroid).  In  contrast,  on  the  broader  system  scale,  beyond 

15  30-km  radius,  the  PV  asymmetry  structure  strongly  resembles  that  in  OBS1  and  suggests 

16  a  left-of-shear  tilt.  Additionally,  an  envelope  of  ascent6  near  30-km  radius  appears 

17  oriented  in  the  right-of-tilt  direction,  as  in  OBS1,  but  is  more  confined  radially.  A 

18  persistent  PV  asymmetry  that  is  approximately  180°  out  of  phase  with  the  near-core  tilt 

19  asymmetry  just  outside  the  peak  heating  radius  (r  =  28  km),  is  likely  generated  and 

20  maintained  by  the  heating  within  this  envelope. 

21 


6  Peak  values  of  diagnosed  ascent  inside  the  region  of  greatest  heating  exceed  40  ms"1  in  OBS5.  The  PV 
perturbation  also  exhibits  local  values  there  that  are  an  order  of  magnitude  larger  than  found  in  OBS1.  In 
the  localized  region  where  YN2— *  0,  these  values  indicate  that  the  limits  of  the  hydrostatic  model  may  have 
been  reached.  Consequently,  the  solution  there  must  be  approached  with  caution. 
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1  5.  Discussion  and  Conclusions 

2  According  to  the  linear  heuristic  model  of  RMG04  for  a  TC-like  vortex  that  is 

3  immersed  in  a  sustained  unidirectional  vertical  shear  flow,  the  rate  at  which  a  vortex 

4  asymptotes  to  the  equilibrium  configuration,  and  hence  its  resilience,  is  a  function  of  the 

5  intrinsic  decay  rate  of  the  core -tilt  mode,  which,  in  the  absence  of  amplification  through 

6  frequency-matched  IGW  radiation,  is  proportional  to  the  mean  radial  PV  gradient  within 

7  the  vortex  skirt.  In  the  context  of  dry  “adiabatic”  linearized  dynamics,  idealized 

8  numerical  solutions  are  shown  to  confirm  that  an  increase  in  the  radial  gradient  of  PV 

9  within  the  vortex  skirt  surrounding  the  core  yields  a  more  rapid  evolution  of  a  sheared 

10  vortex  towards  a  quasi-equilibrium,  left-of-shear  tilt  configuration. 

11  To  further  assess  the  strengths  and  weaknesses  of  the  heuristic  model  in  the  more 

12  realistic  case  of  a  cloudy  vortex  in  unidirectional  shear,  a  simple  parameterization  of 

13  latent  heating  in  the  core  region  is  employed  based  on  prior  work.  For  modest  reductions 

14  of  the  azimuthal-mean  static  stability  within  the  core,  the  heuristic  model  continues  to 

15  provide  a  consistent  picture  of  the  simulated  vortex  evolution  and  its  sensitivities. 

16  However,  in  the  limit  of  strictly  moist-neutral  dynamics  within  a  narrow  region  of  the 

17  “eyewall”,  the  heuristic  model  yields  a  poor  approximation  to  the  simulated  vortex 

18  evolution  within  the  region  of  peak  heating,  but  a  left-of-shear  tilt  of  the  near-core  vortex, 

19  radially  beyond  the  heating  region,  remains  the  preferred  long-time  solution. 

20  The  foregoing  findings  naturally  motivate  the  following  question:  Is  the  heuristic 

2 1  model  still  applicable  to  the  vortex  tilt  evolution  in  the  limit  of  moist-neutral  dynamics, 

22  and,  if  so,  in  what  sense? 


25 


1 


To  help  answer  this  question,  Fig.  12  summarizes  the  transition  in  shear-forced 

2  asymmetry  evolution  from  the  limit  of  no  heating  (OBS1,  Fig.  12a)  to  the  limit  of  zero 

3  moist  stability  (OBS5,  Fig.  12d)  within  the  core.  Consistent  with  our  analysis  in  section 

4  4,  the  PV  perturbation  in  OBS1  (Fig.  12a)  clearly  indicates  the  presence  of  a  forced, 

5  damped  tilt  mode  with  an  amplitude  that  peaks  in  the  core  and  decays  radially  within  the 

6  skirt.  Careful  examination  of  these  solutions  shows  also  radially-outward  propagating 

7  features.  Reasor  and  Montgomery  (2001)  identified  similar  outward-propagating  features 

8  in  the  free-alignment  context  as  trailing  spiral  VRWs.  Since  the  tilt  quasimode  solution 

9  resides  in  the  continuous  spectrum,  it  must  ultimately  decay  through  the  “leaking”  of 

10  energy  via  the  radial  propagation  of  these  VRWs.  They  noted  that  nonlinear  advection 

1 1  tends  to  counteract  the  VRW  dispersion.  When  the  static  stability  is  halved  locally 

12  (OBS2,  Fig.  12b),  the  forced  modal  structure  is  still  evident,  but  the  outward  energy 

13  propagation  is  more  prominent  in  the  solution.  For  small,  but  non-zero,  local  values  of 

14  static  stability  (OBS4,  Fig.  12c),  the  presence  of  a  forced  mode  is  not  readily  discerned  in 

15  the  core.  Instead,  the  PV  perturbation  amplitude  peaks  where  the  heating  is  maximized  (r 

16  ~  28  km).  In  the  limit  of  locally  neutral  (zero)  static  stability  (OBS5,  Fig.  1 2d),  a  similar, 

17  but  larger-amplitude  and  more  radially  confined  PV  perturbation  is  evident  near  the 

18  radius  of  peak  heating.  This  amplitude  maximum,  observed  for  small  values  of  local 

19  static  stability,  is  tied  to  the  vortex  tilt  and  is  likely  generated  and  maintained  by  the 

20  heating  there.  Unique  to  the  moist-neutral  limit,  a  second  and  distinct  region  of  large - 

21  amplitude  PV  asymmetry  exists  radially  inside  the  heating  maximum.  As  fine-scale 

22  banded  features  develop  here  (see  Fig.  11),  Fig.  12d  shows  that  the  PV  perturbation 

23  maintains  its  large  amplitude.  Nonlinear  dynamics  may  become  important  in  this 
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1  “eyewall”  region,  requiring  a  future  modification  of  the  linear  theory.  In  both  OBS4  and 

2  OBS5  a  PV  maximum  outside  30-km  radius  is  associated  with  the  trailing  spiral  VRWs  at 

3  the  outer  edge  of  the  core. 

4  The  timing  and  outward -propagation  speed  of  the  VRWs  are  similar  in  all  OBS 

5  simulations.  Comparison  of  the  simulated  solution  to  the  heuristic  model  solution  in 

6  section  4  reveals  that  a  tilt  mode  continues  to  govern  the  longer-time  evolution  of  vortex 

7  tilt  outside  the  core  in  OBS4  and  likely  in  OBS5.  Thus,  while  the  simulated  linear 

8  primitive  equation  solution  outside  the  heating  region  exhibits  continuous-spectrum  type 

9  VRW  characteristics,  the  tilt  evolution  in  this  part  of  the  near-core  region  is  still 

10  meaningfully  described  by  a  shear- forced  discrete  tilt  mode.  An  assessment  of  whether 

1 1  the  tilt  mode  is  forced  at  all  near  the  core  radius  rc  requires  the  application  of  semi- 

12  analytical  methods,  as  in  Schecter  et  al.  (2002),  which  we  reserve  for  future  work. 

13 
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1  List  of  Figures 

2  FIG.  1.  Illustration  of  the  core  and  skirt  regions  of  an  idealized  tropical  cyclone-like 

3  vortex.  The  region  of  greatest  PV  is  bounded  by  the  core  radius  rc,  and  is  surrounded  by  a 

4  skirt  of  much  smaller  PV.  Here,  rc  is  defined  as  the  radius  of  maximum  PV  gradient.  The 

5  core  may  support  a  discrete  vortex -Rossby  wave  which  resonates  with  the  surrounding 

6  flow  at  the  critical  radius  rt.  Adapted  from  Fig.  2  of  Schecter  et  al.  (2002). 

7  FIG.  2.  The  azimuthal-mean  vortex  used  in  (a)  RWS1  (solid),  RWS2  (dashed)  and  RWS3 

8  (dotted),  and  in  (b)  OBS1-5.  The  inset  in  (a)  shows  the  radial  profile  of  relative  vertical 

9  vorticity  within  the  vortex  skirt,  highlighting  differences  in  the  corresponding  radial 

10  gradient.  For  the  equivalent  barotropic  system,  the  azimuthal-mean  PV  is  related  to  the 

1 1  relative  vorticity  by  q  —  N2(f  +  C, ) . 

12  FIG.  3.  Vortex  alignment  as  a  function  of  skirt  PV  gradient.  Log-linear  plot  shows  the 

13  core  VRW  geopotential  amplitude  evaluated  at  the  core  radius  rc  as  a  function  of  time  for 

14  RWS1  (thick  solid),  RWS2  (dashed)  and  RWS3  (dotted).  The  light  solid  gray  lines 

15  indicate  linear  best  fits  to  the  data  in  each  case.  Time  is  in  units  of  circulation  period  at  rc ; 

16  I0i  I  max  is  the  maximum  value  of  the  geopotential  perturbation  over  each  simulation. 

17  FIG.  4.  Tilt  amplitude  and  direction  in  the  shear-forced  (a)  RWS1,  (b)  RWS2  and  (c) 

18  RWS3  experiments.  Tilt  amplitude  is  measured  by  the  PV  perturbation  (s'3)  at  the  core 

19  radius  rc  (top).  Tilt  direction  is  defined  as  the  crest  azimuth  of  the  PV  perturbation  at  z  = 

20  H  (bottom).  The  solid  lines  show  the  simulated  solution  to  the  primitive  equation  model 

21  of  section  2,  and  the  dashed  lines  show  the  heuristic  model  approximation  to  it  using  the 

22  decay  rate  and  precession  frequency  from  the  respective  free-alignment  simulations.  The 

23  heuristic  model  amplitude  is  scaled  so  as  to  match  the  initial  peak-simulated  value.  90° 
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1  represents  the  left-of-shear  direction,  the  equilibrium-tilt  direction  of  the  sheared  vortex 

2  in  all  cases.  Time  is  in  units  of  circulation  period  at  rc. 

3  FIG.  5.  The  development  of  the  linear  critical  layer  in  RWS3.  Shown  is  the  PV 

4  perturbation  evolution  (from  left  to  right)  within  the  vortex  skirt  for  the  (a)  firee- 

5  alignment  and  (b)  sheared-forced  cases.  The  dashed  line  indicates  the  location  of  the 

6  critical  radius  in  the  free-alignment  experiment.  Here,  qx  max  is  the  maximum  PV 

7  perturbation  over  the  entire  domain  and  over  the  duration  of  each  simulation  (i.e.,  the 

8  peak  core  VRW  amplitude).  Time  is  in  units  of  circulation  period  at  the  core  radius  rc. 

9  FIG.  6.  Temporal  evolution  of  PV  perturbation  amplitude  within  the  vortex  skirt  in  the 

10  shear-forced  RWS3  experiment.  The  dashed  line  indicates  the  location  of  the  critical 

1 1  radius  rt  in  the  corresponding  free-alignment  experiment.  The  modulation  of  disturbance 

12  PV  amplitude  within  the  skirt  occurs  as  a  result  of  the  shear  forcing.  Here,  |qilmax  is  the 

13  maximum  PV  perturbation  over  the  entire  domain  and  over  the  duration  of  the  simulation 

14  (i.e.,  the  peak  core  VRW  amplitude).  Time  is  in  units  of  circulation  period  at  the  core 

1 5  radius  rc. 

16  FIG.  7.  Radial  profdes  of  the  buoyancy  reduction  factor,  Y  for  OBS1  (solid  red),  OBS2 

17  (dotted),  OBS3  (dashed),  OBS4  (dash-dot),  and  OBS5  (solid  black).  For  reference,  the 

18  heavy  solid  line  shows  the  relative  vorticity  profile  of  the  azimuthal-mean  vortex.  In  the 

19  experiments  with  heating,  the  effective  static  stability  is  reduced  at  the  edge  of  the  vortex 

20  core  and  approaches  rapidly  a  typical  “environmental”  value  away  from  the  region  of 

2 1  strongest  tangential  wind. 

22  FIG.  8.  Vortex  alignment  as  a  function  of  effective  static  stability.  Log-linear  plot  shows 

23  the  geopotential  amplitude  as  a  function  of  time  for  OBS 1  (solid  red),  OBS2  (dotted), 
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1  0BS3  (dashed),  0BS4  (dash-dot),  and  OBS5  (solid  black).  In  (a)  the  asymmetry 

2  amplitude  is  measured  at  the  core  radius  rc.  In  (b)  the  asymmetry  is  sampled  at  r  =  40  km, 

3  radially  outside  the  region  of  peak  heating  (cf.  Fig.  7).  The  light  solid  gray  lines  indicate 

4  linear  best  fits  to  the  data  in  each  case.  Time  is  in  units  of  circulation  period  at  rc;  |0ilmax 

5  is  the  maximum  value  of  the  geopotential  perturbation  over  each  simulation. 

6  FIG.  9.  Tilt  amplitude  and  direction  in  the  shear-forced  (a)  OBS1,  (b)  OBS2,  (c)  OBS3, 

7  (d)  OBS4,  and  (e)  OBS5  experiments.  Tilt  amplitude  (top)  is  measured  by  the  PV 

8  perturbation  (s'  )  at  the  core  radius  rc  in  OBS1-3  and  at  r  =  40  km,  radially  outside  the 

9  region  of  peak  heating  (cf.  Fig.  7),  in  OBS4-5.  Tilt  direction  is  defined  as  the  crest 

10  azimuth  of  the  PV  perturbation  at  z  =  H  (bottom).  The  solid  lines  show  the  simulated 

1 1  solution  to  the  primitive  equation  model  of  Section  2,  and  the  dashed  lines  show  the 

12  heuristic  model  approximation  to  it  using  the  decay  rate  and  precession  frequency  from 

13  the  respective  free-alignment  simulations.  The  heuristic  model  amplitude  is  scaled  so  as 

14  to  match  the  initial  peak-simulated  value.  Because  the  precession  frequency  and  decay 

15  rate  cannot  be  meaningfully  estimated  from  the  free-alignment  solution  in  OBS5,  values 

16  from  OBS4  are  used  in  the  heuristic  model  solution  in  (e).  90°  represents  the  left-of-shear 

17  direction,  the  equilibrium-tilt  direction  of  the  sheared  vortex  in  all  cases.  Time  is  in  units 

18  of  circulation  period  at  rc. 

19  FIG.  10.  Time  evolution  of  core  region  PV  perturbation  (shaded,  at  z  =  0)  and  positive 

20  vertical  motion  perturbation  (at  z  =  HI 2)  contoured  at  0.25,  0.5,  0.75,  and  1.5  ms'1  for  the 

21  shear- forced  OBS1  experiment.  Vertical  wind  shear,  indicated  by  the  arrow  in  the  top-left 

22  panel,  is  directed  to  the  right.  The  fields  are  shown  at  (a)  0.25,  (b)  1.5,  (c)  2,  and  (d)  5 

23  circulation  periods.  The  early  downshear-tilted  vortex  corresponds  to  a  PV  dipole  with 
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1  positive  (negative)  perturbation  values  upshear  (downshear)  at  low  levels.  The  horizontal 

2  scale  is  in  km. 

3  FIG.  1 1.  As  in  Fig.  10,  but  for  OBS5.  Peak  values  of  vertical  motion  and  PV  perturbation 

4  within  the  core  exceed  the  contour  ranges  used  in  Fig.  10.  To  facilitate  comparison  of 

5  structure  outside  the  core,  the  contour  levels  used  in  Fig.  10  are  still  employed  here.  See 

6  text  for  further  discussion. 

7  FIG.  12.  Temporal  evolution  of  the  logarithmically-scaled  PV  amplitude  within  the  core 

8  region  in  the  shear- forced  (a)  OBS1,  (b)  OBS2,  (c)  OBS4,  and  (d)  OBS5  experiments. 

9  Time  is  in  units  of  circulation  period  at  the  core  radius  rc. 

10 
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1  Table  1.  Description  of  the  experiments  perfonned  for  this  study.  The  Rankine-with- 

2  skirt  (RWS)  and  observation-based  (OBS)  vortex  profiles  are  shown  in  Fig.  2.  See  text 

3  for  the  definition  of  the  buoyancy  reduction  factor,  Y. 


Experiment 

Vortex  Profile 

Heating  (min  Y) 

RWS1 

RWS:  Skirt  PV  gradient,  control 

No  (1) 

RWS2 

",  t  33% 

No  (1) 

RWS3 

",  t  ioo% 

No  (1) 

OBS1 

Observation-based  PV  monopole 

No  (1) 

OBS2 

66 

Yes  (0.5) 

OBS3 

66 

Yes  (0.25) 

OBS4 

66 

Yes  (0.0625) 

OBS5 

66 

Yes  (0) 
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1  Table  2.  Estimated  vortex  core  radius  rc,  resonant  damping  critical  radius  r„  net  decay 

2  rate  y,  and  core  VRW  precession  frequency  cop  for  each  experiment.  The  net  decay  rate 

3  and  precession  frequency  are  evaluated  at  rc,  except  in  OBS4-5  where  they  are  evaluated 

4  radially  outside  the  region  of  peak  heating  at  r  =  40  km.  The  simulated  net  decay  rate  is 

5  the  sum  of  the  (negative)  resonant  damping  rate  and  (positive)  inertia-gravity  wave 

6  pumping  rate.  See  text  for  further  details.  Because  the  amplitude  decay  in  OBS5  is  not 

7  clearly  exponential  in  time  following  the  initial  adjustment  period  (cf.  Fig.  8),  the  flagged 

8  values  ( ) 1  are  only  tenuously  associated  with  resonant  damping. 


Experiment 

Core  Radius 
rc  (km) 

Critical  Radius 
r,  (km) 

Net  Decay  Rate 
y  (10‘5  s'1) 

Precession 
Frequency 
(Dp  ( 1  O'4  s'1) 

RWS 1/2/3 

30 

69/69/68.5 

-0.10/0.61/1.69 

2.61/2.64/2.70 

OBS 1/2/3 

23 

66.5/64.25/62.75 

-1.97/4.27/6.41 

5.48/5.85/6.13 

OBS4 

23 

61.5 

-10.52 

6.39 

OBS5 

23 

(29)f 

(-5.20)t 

(22.5  l)f 

9 

10 
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1 

2  FIG.  1.  Illustration  of  the  core  and  skirt  regions  of  an  idealized  tropical  cyclone-like 

3  vortex.  The  region  of  greatest  PV  is  bounded  by  the  core  radius  rc,  and  is  surrounded  by  a 

4  skirt  of  much  smaller  PV.  Here,  rc  is  defined  as  the  radius  of  maximum  PV  gradient.  The 

5  core  may  support  a  discrete  vortex -Rossby  wave  which  resonates  with  the  surrounding 

6  flow  at  the  critical  radius  rt.  Adapted  from  Fig.  2  of  Schecter  et  al.  (2002). 
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a)  b) 


x  io"3  x  10 


2  FIG.  2.  The  azimuthal-mean  vortex  used  in  (a)  RWS1  (solid),  RWS2  (dashed)  and  RWS3 

3  (dotted),  and  in  (b)  OBS1-5.  The  inset  in  (a)  shows  the  radial  profde  of  relative  vertical 

4  vorticity  within  the  vortex  skirt,  highlighting  differences  in  the  corresponding  radial 

5  gradient.  For  the  equivalent  barotropic  system,  the  azimuthal-mean  PV  is  related  to  the 

6  relative  vorticity  by  q  =  N2[f  +  C, ) . 
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2  FIG.  3.  Vortex  alignment  as  a  function  of  skirt  PV  gradient.  Log-linear  plot  shows  the 

3  core  VRW  geopotential  amplitude  evaluated  at  the  core  radius  rc  as  a  function  of  time  for 

4  RWS1  (thick  solid),  RWS2  (dashed)  and  RWS3  (dotted).  The  light  solid  gray  lines 

5  indicate  linear  best  fits  to  the  data  in  each  case.  Time  is  in  units  of  circulation  period  at  rc; 

6  I0i  I  max  is  the  maximum  value  of  the  geopotential  perturbation  over  each  simulation. 
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FIG.  4.  Tilt  amplitude  and  direction  in  the  shear-forced  (a)  RWS1,  (b)  RWS2  and  (c) 
RWS3  experiments.  Tilt  amplitude  is  measured  by  the  PV  perturbation  (s'  )  at  the  core 
radius  rc  (top).  Tilt  direction  is  defined  as  the  crest  azimuth  of  the  PV  perturbation  at  z  = 
H  (bottom).  The  solid  lines  show  the  simulated  solution  to  the  primitive  equation  model 
of  section  2,  and  the  dashed  lines  show  the  heuristic  model  approximation  to  it  using  the 
decay  rate  and  precession  frequency  from  the  respective  free-alignment  simulations.  The 
heuristic  model  amplitude  is  scaled  so  as  to  match  the  initial  peak-simulated  value.  90° 
represents  the  left-of-shear  direction,  the  equilibrium-tilt  direction  of  the  sheared  vortex 
in  all  cases.  Time  is  in  units  of  circulation  period  at  rc. 
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FIG.  5.  The  development  of  the  linear  critical  layer  in  RWS3.  Shown  is  the  PV 
perturbation  evolution  (from  left  to  right)  within  the  vortex  skirt  for  the  (a)  firee- 
alignment  and  (b)  sheared-forced  cases.  The  dashed  line  indicates  the  location  of  the 
critical  radius  r*  in  the  free-alignment  experiment.  Here,  qlm ax  is  the  maximum  PV 

perturbation  over  the  entire  domain  and  over  the  duration  of  each  simulation  (i.e.,  the 
peak  core  VRW  amplitude).  Time  is  in  units  of  circulation  period  at  the  core  radius  rc. 
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FIG.  6.  Temporal  evolution  of  PV  perturbation  amplitude  within  the  vortex  skirt  in  the 
shear- forced  RWS3  experiment.  The  dashed  line  indicates  the  location  of  the  critical 
radius  i\  in  the  corresponding  free-alignment  experiment.  The  modulation  of  disturbance 

PV  amplitude  within  the  skirt  occurs  as  a  result  of  the  shear  forcing.  Here,  |<?ilmax  is  the 
maximum  PV  perturbation  over  the  entire  domain  and  over  the  duration  of  the  simulation 
(i.e.,  the  peak  core  VRW  amplitude).  Time  is  in  units  of  circulation  period  at  the  core 
radius  rc. 
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FIG.  7.  Radial  profiles  of  the  buoyancy  reduction  factor,  Y  for  OBS1  (solid  red),  OBS2 
(dotted),  OBS3  (dashed),  OBS4  (dash-dot),  and  OBS5  (solid  black).  For  reference,  the 
heavy  solid  line  shows  the  relative  vorticity  profile  of  the  azimuthal-mean  vortex.  In  the 
experiments  with  heating,  the  effective  static  stability  is  reduced  at  the  edge  of  the  vortex 
core  and  approaches  rapidly  a  typical  “environmental”  value  away  from  the  region  of 
strongest  tangential  wind. 
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FIG.  8.  Vortex  alignment  as  a  function  of  effective  static  stability.  Log-linear  plot  shows 
the  geopotential  amplitude  as  a  function  of  time  for  OBS1  (solid  red),  OBS2  (dotted), 
OBS3  (dashed),  OBS4  (dash-dot),  and  OBS5  (solid  black).  In  (a)  the  asymmetry 
amplitude  is  measured  at  the  core  radius  rc.  In  (b)  the  asymmetry  is  sampled  at  r  =  40  km, 
radially  outside  the  region  of  peak  heating  (cf.  Fig.  7).  The  light  solid  gray  lines  indicate 
linear  best  fits  to  the  data  in  each  case.  Time  is  in  units  of  circulation  period  at  rc ;  |0ilmax 
is  the  maximum  value  of  the  geopotential  perturbation  over  each  simulation. 
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FIG.  9.  Tilt  amplitude  and  direction  in  the  shear-forced  (a)  OBS1,  (b)  OBS2,  (c)  OBS3, 
(d)  OBS4,  and  (e)  OBS5  experiments.  Tilt  amplitude  (top)  is  measured  by  the  PV 

•5 

perturbation  (s'  )  at  the  core  radius  rc  in  OBS1-3  and  at  r  =  40  km,  radially  outside  the 
region  of  peak  heating  (cf.  Fig.  7),  in  OBS4-5.  Tilt  direction  is  defined  as  the  crest 
azimuth  of  the  PV  perturbation  at  z  =  H  (bottom).  The  solid  lines  show  the  simulated 
solution  to  the  primitive  equation  model  of  Section  2,  and  the  dashed  lines  show  the 
heuristic  model  approximation  to  it  using  the  decay  rate  and  precession  frequency  from 
the  respective  free-alignment  simulations.  The  heuristic  model  amplitude  is  scaled  so  as 
to  match  the  initial  peak-simulated  value.  Because  the  precession  frequency  and  decay 
rate  cannot  be  meaningfully  estimated  from  the  free-alignment  solution  in  OBS5,  values 
from  OBS4  are  used  in  the  heuristic  model  solution  in  (e).  90°  represents  the  left-of-shear 
direction,  the  equilibrium-tilt  direction  of  the  sheared  vortex  in  all  cases.  Time  is  in  units 
of  circulation  period  at  rc. 
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FIG.  9.  (cont’d) 
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FIG.  10.  Time  evolution  of  core  region  PV  perturbation  (shaded,  at  z  =  0)  and  positive 
vertical  motion  perturbation  (at  z  =  HI 2)  contoured  at  0.25,  0.5,  0.75,  and  1.5  ms'1  for  the 
shear-forced  OBS1  experiment.  Vertical  wind  shear,  indicated  by  the  arrow  in  the  top-left 
panel,  is  directed  to  the  right.  The  fields  are  shown  at  (a)  0.25,  (b)  1.5,  (c)  2,  and  (d)  5 
circulation  periods.  The  early  downshear-tilted  vortex  corresponds  to  a  PV  dipole  with 
positive  (negative)  perturbation  values  upshear  (downshear)  at  low  levels.  The  horizontal 
scale  is  in  km. 
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FIG.  1 1.  As  in  Fig.  10,  but  for  OBS5.  Peak  values  of  vertical  motion  and  PV  perturbation 
within  the  core  exceed  the  contour  ranges  used  in  Fig.  10.  To  facilitate  comparison  of 
structure  outside  the  core,  the  contour  levels  used  in  Fig.  10  are  still  employed  here.  See 
text  for  further  discussion. 
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2  FIG.  12.  Temporal  evolution  of  the  logarithmically-scaled  PV  amplitude  within  the  core 


3  region  in  the  shear-forced  (a)  OBS1,  (b)  OBS2,  (c)  OBS4,  and  (d)  OBS5  experiments. 


4  Time  is  in  units  of  circulation  period  at  the  core  radius  rc. 
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